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Abstract 

The properties of the low excitation field magnetic response of the granular high temperature (HTc) superconductor 
Lai,8Sro.2Cu04 have been analyzed at low temperatures. The response of the Josephson currents has been extracted 
from the data. It is shown that intergrain current response is fully irreversible, producing shielding response, but 
do not carry Meissner magnetization. Analysis of the data shows that the system of Josephson currents freezes into 
a glassy state even in the absense of external magnetic field, which is argued to be a consequence of the d-wave 
nature of superconductivity in Lai.8Sro.2Cu04. The macroscopic diamagnetic response to very weak variations of the 
magnetic field is shown to be strongly irreversible but still qualitatively different from any previously known kind of 
the critical-state behaviour in superconductors. A phenomenological description of these data is given in terms of a 
newly proposed "fractal" model of irreversibility in superconductors. 
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I Introduction 

Granular superconductors (SC) are composed of a very 
large number of small (micron-size) superconductive 
grains which are coupled together due to the Joseph- 
son tunnelling (or, in some cases, due to the proxim- 
ity effect). These systems are inherently disordered due 
to randomness in the sizes of grains and in their mu- 
tual distances. Usually the strength of Josephson cou- 
pling between grains is rather weak, so the maximum 
Josephson energy of the contact between two grains is 
much below the intragrain superconductive condensa- 
tion energy. Therefore granular SC can be considered 
as systems with a two-level organization; their short- 
scale properties are determined by the superconductiv- 
ity of individual grains, whereas the macroscopic SC be- 
haviour is governed by the weak intergrain couplings. In 
the treatment of the latter, one can neglect any inter- 
nal structure of SC grains and describe them just by 
the phases of their superconductive order parameters 
Aj = \A\j exp{i(j>j). As a result, the macroscopic be- 
haviour of granular SC can be described by a classical 
free energy functional of the form (cf. Ref. 0, ||, |3|): 



(fr i — (rotAf - —rot A ■ Ue^t 



(I.l) 



where a^- = ^ /-"' Adr is the phase difference induced by 
the electromagnetic vector potential A and $o = nhc/e, 
whereas the coupling strengths Ej are proportional to 



the maximum Josephson currents: Sy — ^Ifj 



The 



vector potential A in Eg. (p|) is the sum of the vector 
potential Aext of the external magnetic field Hext and of 
the Josephson currents- induced vector potential Aind- 
In the absence of external magnetic field, the lowest- 
energy state for the "Hamiltonian" (^t|) is, clearly, a 
macroscopically superconductive state with all phases 
equal to each other. Thus that granular SC system looks 
similar to the random XY ferromagnet with randomness 
in the values of the coupling strengths Ej 's (apart from 
the possible role of the induced vector potential A^d 
which will be dicussed later); within this analogy the 
role of XY "spin components" is taken by Sx = cos , 
Sy = sin . 

The situation becomes a lot more complicated in 
the presence of non-zero external magnetic field Hext, 
which makes the system randomly frustrated (since mag- 
netic fluxes penetrating plaquettes between neighbouring 
grains are random fractional parts of $o)- When the ex- 
ternal field is sufficiently strong, Hext ^ Hq = $o/«o 



(here oq is the characteristic intergrain distance), the 
random phases aij become of the order of n or larger, 
which means complete frustration of the intergrain cou- 
phngs - i.e. the system is then expected to resemble the 
XY spin glass. Actually the random Josephson network 
in a magnetic field is not exactly identical to the XY 
spin glass due to the following reasons §: i). The effec- 
tive couplings Ej = Ej exp(iaij) between "spins" Si of 
the frustrated SC network are random complex numbers 
whereas in the XY spin glass model, they are real random 
numbers, ii). Generally the phases aij depend on the to- 
tal magnetic induction B — Hext+'^ind^ i-e. the effective 
couplings Ej depend on the phase variables (j)j determin- 
ing the intergrain currents = IfjSin{(j)i — 4>j — aij). 
In some cases the effects produced by the self-induced 
magnetic field 'Bind are weak and can be neglected (the 
quantitative criterion will be discussed later on) , so that 
phases aij can be considered as being fixed by the exter- 
nal field. 



The model described by the Hamiltonian (I.l) with 



fixed aij 's and Hext ^ Hq is usually called "gauge glass" 
model. It is expected on the basis of the analytical |§, ^, 
^, H as well as numerical ^ results that the gauge glass 
model in 3D space exhibits a true phase transition into 
a low-temperature glassy superconductive (nonergodic) 
state. The mean-field theory of such a low-temperature 
state shows || that it is characterized by the presense 
of a finite effective penetration depth for the variation of 
an external field, nonzero macroscopic critical current, 
and the absense of a macroscopic Meissner effect. The 



full model ([.1) with a's containing contribution from 
Bind is sometimes called "gauge glass with screening" Q . 
The effect of screening on the presence and properties of 
the phase transition into a glassy state is not completely 
clear; some numerical results |^ indicate the absence of 
a true phase transtion in a 3D model with screening. 
Quantitatively, the strengh of screening is determined by 
the ratio Pl = 27r£/c/c<I>o where C is the characteristic 
inductance of an elementary intergrain current loop 
In the ceramics with /J^ ^ 1 , screening effects become 
important on a long-distance scale ~ ciq/^/]3l only (i.e. 
they are similar to the strongly type-II superconductors 
with disorder). 

Apart from its relevance for the description of granular 
superconductors, the gauge glass model with screening 
is rather often considered (e.g. Ref. ||l^) as a simpli- 
fied model describing the large-scale behaviour of disor- 
dered bulk type-II superconductors in the mixed state 
(so-called vortex glass problem). Actually it is unclear 
a priori how these two problems are related; an obvi- 
ous difference between them is that the basic ingredi- 
ent of the latter is the vortex lattice which is clearly 
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an anisotropic object, whereas the former does not con- 
tain any prescribed direction in the 3D space. On the 
other hand, the granular superconductor in a moderate 
magnetic field H^xt < niay be considered as a kind 
of disordered type-II superconductor, where the notion 
of a hypervortex (which is the macroscopic analogue of 
the Abrikosov vortex) can be introduced ^, |ll| . There- 
fore, the macroscopic properties of a granular network 
at Hext ^ ^^0 may resemble those of the vortex glass; in 
such a scenario a phase transtion between vortex glass 
and gauge glass phases would be expected in a granular 
superconductive network at Hext ~ Hq (cf. Ref.[|l| for a 
more detailed discussion of this subject). 

Recently, it was noted that granular superconductors 
may become glassy even in the absence of external mag- 
netic field, if a large enough part of Josephson junctions 
are anomalous, i.e. their minimum Josephson coupling 
energy corresponds to a phase difference Acj) — tt instead 
of (so-called 7r-j unctions). Two completely different 
origins of 7r-j unctions were proposed: mesoscopic fluctu- 
ations in dirty superconductors 1 12 and the pairing with 
non-zero momentum |l3| , [l^ . Recent experiments re- 
vealing the d-wave nature of pairing in high-temperature 
superconductors (l5) indicate the possibility of observing 
glassy superconductive behaviour in HTSC ceramics in 
virtually zero magnetic field. Note that ceramics with 
equal concentrations of usual and 7r-junctions are com- 
pletely equivalent (if screening effects can be neglected) 
to the XY spin glass. Contrary to the 3D gauge glass 
model, the XY spin glass in 3D is expected to have no 
true thermodynamic phase transition at finite tempera- 
ture 1^ ; recently, it has been suggested that the XY spin- 
glass and d-wave ceramic superconductor might have a 
new equilibrium ordered phase, the so-called chiral-glass 
phase [|l6j. However, these issues are hardly relevant for 
the measurable response at temperatures much below 
"bare" glass transition temperature Tg, which we con- 
sider in this paper. 

Experimental studies of granular superconductors re- 
veal an appearance of magnetic irreversibility (a 
difference between Meissner and shielding magnetiza- 
tions or, in other terms, between Field Cooled (F.C.) and 
Zero Field Cooled (Z.F.C.) magnetizations) below some 
temperature Tg, which is lower than the SC transition 
temperature Tc of the grains. However, detailed analysis 
of the magnetic response in such systems is usually com- 
plicated by the mixing of contributions from individual 
grains and from the intergrain currents. The goal of this 
paper is to develop a method which makes it possible 
to extract from the raw data on d.c. magnetic response 
the intergrain contribution and to compare its behaviour 
with existing theoretical predictions. 



The compound Lai.8Sro.2Cu04 was chosen in this 
study for experimental convenience, since its critical tem- 
perature (« 32 K ) is within the optimal temperature 
range of our noise and a.c. susceptibility measurements 
setup. The sample was fabricated by standard solid state 
reaction of La203, SrCOs and CuO Mixed powder 
was pressed into pellets which were sintered in air at 
920°C for 12 hours. The material was then submitted to 
three cycles of regrinding, sifting to 20 fj,m, pressing and 
sintering again at 1100°C for 12 hours. Samples prepared 
in two successive runs were used in this study. In the first 
one (sample A), pellets 1 mm thick and 10mm diameter 
were obtained, with a density about 80% of the theoreti- 
cal bulk value. In the second one (sample B), cylinders of 
diameter 6 mm and length 5 to 6 mm were prepared with 
a density ratio about 88%. In both preparations, grains 
sizes were in the range 1-10 /im. Room temperature 
X-ray powder diffraction patterns showed the presence 
of a small amount (< 5%) of the non superconductive 
compound Lai^2a;Sr2a;Cu205. 

The rest of the paper is organized as follows. In Sec- 
tion II the general analysis of the magnetic response 
data obtained on two different samples (A and B) of 
Lai.8Sro.2Cu04 ceramics is presented and the intergrain 
(Josephson) contribution to the overall response is ex- 
tracted. Section HI is devoted to the detailed study of 
the magnetic response of Josephson intergrain network in 
the low-field range. It is found that the macroscopic crit- 
ical current is suppressed considerably (by a factor 2), in 
a magnetic field of only about 2G. The lower-field d.c- 
response to field variations of order 0.05 0.5G was ana- 
lyzed for the F.C. states obtained at Hpc — 0-f- lOG and 
two temperatures, lOK and 20K. The data at T = lOK 
and Hpc — and O.IG are shown to be compatible 
with the Bean critical-state picture ||l9| and the low- 
field critical current value is identified. The rest of the 
data are in a sharp contrast with Bean-model predic- 
tions: the screening current grows suhlinearly (approxi- 
mately as a square root) with increasing exitation field. 
Very low field, low frequency a.c. measurements are 
presented, which reveal the strongly irreversible nature 
of that anomalous response. A new phenomenological 
model is proposed for the treatment of these data. Its 
first predictions are found to be in a reasonable agree- 
ment with the data. In Section IV, the theoretical analy- 
sis of our experimental results is given in terms of the ex- 
isting theories of "gauge-glass" state. It is shown that the 
observed transition temperature to the low-temperature 
state of the network and the magnitude of the (low- 
B, T) critical current are in sharp contradiction with the 
(usual) assumption that the zero-field granular network 
is unfrustrated. On the contrary, under the assump- 
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tion of a strongly frustrated network a,t B — 0, all basic 
measured parameters of the ceramic network are in mu- 
tual agreement. We believe that these estimates indicate 
the existence of a large proportion of 7r-junctions in the 
Lai.8Sro.2Cu04 ceramics, possibly due to the d-wave 
nature of superconductivity in cuprates. The Section V 
is devoted to the development of a new model of dia- 
magnetic response in glassy superconductors, which is 
necessary for the description of the anomalous data de- 
scribed at the end of Section III. This new model (in 
some sense, intermediate between the Bean and the 
Campbell ones) is based on two ideas: i) the ex- 
istence of two characteristic "critical" currents ( Jd and 
Jc ^ Jci), and ii) the fractal nature of free energy valleys 
in the ceramic network. Our conclusions are presented 
in Section VI, whereas some technical calculations can 
be found in the Appendix. 

For convenience, the e.m.u system of units will be used 
for experimental data, and Gaussian units for the theo- 
retical discussions. 



the field parallel to the longitudinal axis is ~ 0.06 |22 



II General Properties 
Magnetic Response 



of D.C. 



The d.c. magnetization was measured by the classical 
extraction method. Two SQUID magnetometers were 
used: one a home made apparatus used in several previ- 
ous spin-glass studies , the other a commercial system 
(Cryogenics S500). 

In this section, we describe successively the static mag- 
netic response of samples A and B and present a prelim- 
inary treatment of these data, in order to distinguish 
between the magnetic response of individual grains and 
intergrain currents (|, 22| (a detailed study of the latter 
is the subject of the next section). Firstly, we present 
results obtained after cooling the samples in various d.c. 
fields and applying small field increases. Secondly, we 
will derive from the results the response of the Josephson 
currents as a function of field and temperature. Finally, 
we will show that the behavior of the field cooled (F.C.) 
susceptibility can be satisfactorily accounted for if the 
system of Josephson currents does not carry Meissner 
magnetization. It will be shown that the same interpre- 
tation accounts fairly well for the F.C. results which, at 
first sight, are rather different for the samples A and B. 

II-A Sample A 

Sample A is a 1 mm thick pellet with an approximately 
ellipsoidal shape of 2 x 6 mm. Its calculated volume is 
V « 8.5 mm^ and the demagnetizing field coefhcient for 



Figure 1: Magnetic moment of the sample A as a func- 
tion of field applied in the zero field cooled state (e.m.u. 
units of moment correspond to cm'^ — G) . 

Fig. |l| displays the magnetic dipole moment of the sam- 
ple cooled to 10 K in zero field and submitted to cycles 
H„iax — > for several values of Hmax up to 2 G. 
At the lowest increasing fields, the moment increases 
initially with a slope Pi. Above 1.5 G, it approaches 
a slope P2. The remanent positive moment saturates 
for Hmax > 1 G. The calculated moment of the sample 
for perfect volume shielding in an homogeneous field is 
(e.m.u. system): 



7W = - 



H-V 
47r(l - N) 



= -0.72- 10"-^ X Hcm^-G. 



Owing to the error in the evaluation of the volume, this 
value is determined with an accuracy of only ±5% . Nev- 
ertheless, it is in fair agreement with the slope PI in 
Fig. |l|. On the other hand, the slope P2 is about 53%, a 
rather small value since the density ratio of the sample 
is about 80%. At such low temperatures (in comparison 
with Tc w 32 K) , where the lower critical field of the 
grain's material is above 100 G, one would expect expul- 
sion of the field by the grains with a penetration depth A. 
The expected value for the magnetization M — M./V oi 
the system of uncoupled grains system can be calculated 
as [g: 

E = / 
H Att' l~ fN - {\- f)n 



{11.2) 



where / is the volume fraction of the superconduc- 
tive material and n is the demagnetizing field for the 
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grains. For an estimate, we assume grains to be spherical 
(n = 1/3) and, using M/H « 0.53 • l/47r and N = 0.06, 
we find / « 0.41. This value is considerably below the 
volume fraction of the sample filled by grains (w 0.8); 
we assume that the difference is due to the intergrain 
penetration depth A being comparable to the grain size 
r and estimate an effective value of A as 

/ ^ 0.41 = 0.8 ^1 - ^ yielding A = 0.2 r . 

Taking an average size of 5 for the grains, we obtain 
A « 500 nm. Values reported for the mean penetration 
depth in Lai.8Sro.2Cu04 are about 200 25|. The 
value found here is larger than the expected mean value 
for the homogeneous material, indicating that the grains 
are not monocrystalline. This will be confirmed below 
by the results of field-cooling experiments. 



Figure 2: Shielding susceptibiliy of the sample A as a 
function of temperature, normalized to the moment for 
complete shielding. Curves are arranged in the same 
ascending order as in the legend. 

The shielding susceptibility is plotted in Fig. |^, as a 
function of temperature and for several values of the 
ambient F.C. field. The measurements were performed 
according to the following procedure: the sample was 
cooled in a field iJ^c down to the working tempera- 
ture and the moment was measured after waiting 300 
sec; then the field was increased by a small amount 
A_ff < Hdc/W and the moment was measured again af- 
ter waiting 300 sec. The figure displays the experimental 
shielding susceptibility normalized to the value for total 
shielding, i.e.: 



The curves show the double step usually ascribed to 
the action of both intragrain currents and Josephson in- 
tergrain currents p6[ . At high temperature, the onset 
of grains diamagnetism occurs at about 32 K. Above 
25 K, the response corresponds to the diamagnetism of 
the grains. At a fixed temperature, it is Hdc inde- 
pendent for Hdc < 5 G, and decreases for increasing 
Hdc > 5 G. Below 25 K, the onset of Josephson currents 
manifests as a second step of the diamagnetic response. 
This second step appears at a decreasing temperature 
as Hdc increases. At the lowest temperatures, the dia- 
magnetic moment amounts to about 100% of flux expul- 
sion at Hdc = and decreases with increasing Hdc- At 
Hdc > 8 G, the flux expulsion saturates at a value slightly 
above 50% which corresponds roughly to the level of 53% 
determined above for the grains response. 

The susceptibility in Fig. ^ contains the contributions 
of grains and Josephson currents. The contributions 
can be separated on the line of the work by Dersh and 
Blatter |Q. The induction in the sample is given by 
B ^ H + 47r(Afg -I- Mj) where Mg and Mj stand re- 
spectively for the magnetization of grains and of the 
Josephson currents. It should be noted that the mag- 
netization due to macroscopic circulating currents in a 
superconductor is sample-size dependent, i.e. the corre- 
sponding susceptibility is not a local quantity. At the 
macroscopic scale of the circulating currents, the mag- 
netization Mg can always be written as XgHiocai, where 
Xg{H) is homogeneous over the sample. In what fol- 
lows, we consider quantities averaged over the volume 
of the sample: in that case, Mj is the averaged mo- 
ment per volume unit due to the currents. The de- 
magnetizing field effect will be neglected in the calcu- 
lations. We have verified that, owing to the small value 
of the demagnetizing factor , this does not modify the 
essential features of the result while allowing a sim- 
pler derivation (the effect of demagnetizing factor will 
be taken into account when analyzing the data from 
the sample B). We get Mg = XgiH + iirMj) . Then 
M = Mg + Mj = XgH + MjUg (with ^lg = 1 + 47rxg ), 
and 

M - XgH 



(II.3) 



Xsh 



M{H + AH) - M{H) ^ 47r(l - N) 



Eq. ( pi.3| ) must be considered with care since Xg is his- 
tory and field dependent. In fact it is well-adapted to 
the description of the result of zero (or small) field cool- 
ing experiments. More generally, we must consider the 
response to field increments 5H to obtain x = 5M/5H . 
Then, the polarizability ||2^ Xj of the Josephson network 
reads: 

X-Xg 



AH 



V 



Xj 



(11.4) 
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Note that we can equivalently consider the response of 
the currents system in an homogeneous medium with 
permeabihty fj,g. If the apphed field is varied by SH, 
the Josephson network sees a variation of internal field 
SHi = iigSH and develops a polarization 5Mj ~ Xj^Hi . 
Then, we recover Eq. (II.4). 

The value of Xg could be determined in principle if 
we were able to obtain a packing of disconnected grains 
equivalent to the packing of the sintered sample. In prac- 
tice this was not possible. Indeed, mechanical grinding 
resulted in breaking a large part of the grains and thus 
modifying the characteristics of the material. Neverthe- 
less, it is possible to extract Xg, at least approximately, 
from the data of Fig. |^. At high temperature, above 
the onset of intergrain currents at « 25 K, the shielding 
susceptibility Xsh is due to the grains alone, indepen- 
dent on Hdc below « 6 C? . At low temperatures, for Hdc 
above « 6 G , the Xsh curves superpose and there is no 
manifestation of the onset of intergrain currents. Thus , 
here also, Xsh represents the response of the grains alone. 
Hence, the response Xg of the grains can be reasonnably 
approximated by an interpolation between these two lim- 
its. The interpolation curve, obtained by a smooth- 
ing procedure between both curves at Hdc — OG and 
Hdc = 20G is displayed on Fig . p| ( dotted curve). The 
values of Xj derived from Eq. (IL4) are plotted versus 
temperature in Fig. 0, for Hdc < 6 G. 



sponse due to the correlative increase of the value of AH 
{AH = H/10) cannot be excluded: this aspect will be 
studied in detail in sample B. Finally, one can note the 
similarity of our data with the results of earlier numeric 
simulations on a gauge glass system . 

Above we have dicussed the system's responses to 
the variation of magnetic field at fixed temperature (i.e. 
shielding responses) and extracted from these data the 
polarizability Xj of the intergrain system. Now we turn 
to the description of the results of the Field Cooling 
(F.C.) measurements. E.G. (Meissner) magnetization 
was measured by the standard procedure between lOK 
and 40 K for fields from 0.01 to 20 G. The results are 
reported in Fig. ^ versus temperature and Fig. |^ versus 
applied field. Data are normalized to the value of the 
moment for 100% shielding. 



Figure 4: F.C. (Meissner) susceptibility of sample A as a 
function of temperature for fields up to 20 G . Curves are 
arranged in the same ascending order as in the legend. 



Figure 3: Josephson currents susceptibility in the sam- 
ple A as a function of temperature. Values have been 
calculated from data of Fig. and using Eq. (II.4). 



Note that the dependence of Xj on Hdc , seen in the 
figure is supposed to reflect the behavior of the initial 
shielding properties of the Josephson network with the 
increase of Hdc- Nevertheless, non linearity of the re- 



Even at the smallest field, the flux expulsion rate is 
no more than 45%, less than the 53% shielding by the 
grains. At low fields, below 1 G, there is an approximate 
affinity between the curves of M/H versus T. M/H can 
be extrapolated linearly to H — > 0. The result is plotted 
in Fig. H (solid circles): one can see that the extrapo- 
lated F.C. susceptibility superposes exactly with the low 
d.c. field shielding susceptibility above 25 K. Therefore, 
at low d.c. field above 25 K, the response of the grains 
system is reversible and it is well described by the low 
d.c. field shielding curves; this justifies the hypothesis 
used above for the calculation of Xj ■ On the other hand 
(see Fig. ||), the behavior of the F.C. susceptibility as 
a function of H is not trivial. M/H decreases with in- 
creasing field and reaches a stable level (about 25% at 
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Figure 5: F.C. (Meissner) susceptibility of sample A as 
a function of field for selected temperatures Curves are 
arranged in the same ascending order as in tthe legend. 



the lowest temperatures) at roughly 1 G. Whatever the 
temperature, this decrease is centered at a constant value 
of the field, about 0.1 — O.SG. Above 5G , M/H decreases 
once more with increasing field. Note an essential differ- 
ence between the F.C. results presented on Fig. ^ and 
the shielding results above (Fig. |^): the F.C. curves do 
not show any increase of the response M/H with the 
temperature decrease below 20K, where the intergrain 
coupling grows considerably (as it is seen from Fig. |2|). 
This means that the network of intergrain currents does 
not produce Meissner (F.C.) magnetization, whereas it 
does produce shielding magnetization. 

The behavior of the F.C. susceptibility xfc — 
Mpc/H as a function of the applied field H depicted 
in Fig. ^ shows two nontrivial features: i) crossover be- 
tween two plateaus (at low and moderate values of H), 
which takes place between O.IG and IG independent of 
temperature, and ii) the value of the low-field xfc is no- 
ticeably lower than the Meissner response of uncoupled 
grains (45% versus 53%). These features can be under- 
stood in terms of (i), a polycrystalline structure of the 
grains, which can be suspected from the large values of 
the penetration depth obtained from the results of Fig. [|, 
and (ii), self-shielding (pinning of the magnetic flux) by 
the Josephson currents when lowering the temperature 
in an applied field. 

We start from the feature i); the curves of F.C. magne- 
tization in Fig. 1^ are rather similar to those which were 
measured by Ruppel et al. |^ in YBaCuO ceramics. 
The authors interpreted their results on the basis of a 
theory of the flux expulsion by strongly anisotopic ran- 



domly oriented crystallites as derived by Wohllebeen et 
al. We stress that the model is not based on any 

activated flux creep mechanism. It is thus well-adapted 
to the analysis of our results: indeed, flux creep effects 
can hardly be invoked here since the temperature has no 
apparent effect of on the characteristic field related to 
the decrease of magnetization. The starting point of the 
model is that, provided the size b of the crystallites is 
such that A|| <C & <IC Aj_, the longitudinal magnetization 
of a crystallite whose c-axis makes an angle a with the 
field is given by M — —{H/Att) ■ 7(cosa)^, where 7 is a 
factor close to 1, depending on the ratio A||/6. After aver- 
aging over a, one obtains M/H = (7/3)(l/47r). It must 
be stressed that the system of intragrain crystallites is a 
strongly-coupled system, contrary to the system of grains 
which composes the ceramic. Therefore, a grain consists 
of an ensemble of interconnected Josephson loops sur- 
rounding crystallites whose planes are nearly along the 
field and are thus transparent to the field. At low fields, 
this system will expel the flux with a penetration depth 
depending on the junction coupling energy. Nevertheless, 
when the field is such that a loop sees a flux larger than 
^ $0/2, the macroscopic magnetization of the Josephson 
currents vanishes and the system reacts as an ensemble 
of disconnected crystallites BOll. The characteristic field 
of this crossover is such that 



H„ 



0.1 



(II.5) 



Recently determined values for the penetration depth in 
Lai gSro.2Cu04 |3l| are Ay = 150 nm and A_l = 1500 nm. 
Older measurements indicate a higher anisotropy, up to a 
factor 14 . We can thus reasonnably consider that the 
model can be applied in our case. Taking Hm = 0.3 G, 
we obtain Sc — 7.4 • 10~^cto^ . With Sc ~ 7r6^ this 
leads to a mean diameter 6 = 1.5 /zm for the crystal- 
lites. Above Hjrn the system acts as an ensemble of 
crystallites whose average susceptibility is (7/3)(l/47r) . 
With the density ratio / = 0.8, taking 7 = 1 and sup- 
posing spherical crystallites we obtain from Eq. ( II. 2 ) 



47r^ = 0.31 which is above the experimental value ( the 
latter being about 0.25). Nevertheless, it must be noted 
that we have neglected here the effect of the factor 7 and 
used a rather unrealistic spherical approximation for the 
shape of crystallites. Finally, it has been seen that above 
5 G , the F.C. magnetization starts to decrease once more 
with increasing field although Hci is larger than lOOG in 
Lai.8Sro.2Cu04 . This can be due to intrinsic pinning 
inside the crystallites themselves when the applied field 
is such that the flux in the cross section of one crystallite 
is larger than <I>o . With a mean radius of 0.8 (xm for the 
crystallites, this crossover occurs at about 10 G . 



6 



Now we turn to the discussion of the feature ii) men- 
tioned above. At temperatures below 25 K , the Joseph- 
son currents become active. Their effect is that, at 10 ii' , 
the shielding response of the system of grains amounts 
at about 53%, while the F.C. susceptibility saturates at 
about 45%. This difference is enough to be significant 
and can be interpreted as the result of pinning by the 
Josephson network. In fact, this pinning can be under- 
stood as a back shielding effect of the Josephson cur- 
rents against the decrease of local internal field, due to 
the temperature dependence of the grain's system per- 
meability /ig. We have seen above that the response of 
the system consists of the two parts: (i) for an applied 
field H, the internal field due to the grains seen by the 
intergrain currents is Hi — fig H , and (ii) the intergrain 
currents system reacts to all variation of Hi with a polar- 
izability Xj and generates a magnetization 6Mj — Xj^Hi. 
Thus, when the temperature is decreased by dT, the in- 
ternal field decreases by Hd^g/dT and the Josephson 
network tends to screen this variation. Since the inter- 
grain currents give no Meissner effect, we consider their 
response as totally irreversible. Thus for a variation dT 
of the temperature, in a field H , the variation of induc- 
tion is: 



dB 



{l+^TTX,)^ 



On the other hand, B = {1 
1 + 47rXg , we finally obtain: 

Xfc = Xg + 47r 



P-) HdT 



With 



Ms 



Xj ^dT^Xg+X^ 



FC 
j 



(11.6) 



Mr 



,FC 



Xj H is the magnetization produced by the 
Josephson currents due to variation of fig with decreasing 
temperature. As Xg is known only in the limit Hdc , 
Eq. (II.6) has been used to calculate xfc versus T in 
the limit of low field. In order to do it, we started from 
the values of Xgi^^ ~^ 0) a^s derived above; for Xj: we 
have used the values given in Fig. || for the smallest field 
Hdc = 0.06 G. The result is plotted on Fig. | The 
agreement of calculated values with experimental data is 
rather satisfactory, although not perfect. 



This discrepancy is emphasized if Eq. ([1.6) is reversed 
in order to calculate Xj as a function of Xg and xfc ■ The 
reason is that we have used here the simplest linear model 
of back shielding. In fact, as we will see later, the re- 
sponse of the currents system is strongly non-linear, with 
the susceptibility Xj decreasing with increasing AH , and 
this effect becomes stronger as the temperature increases. 
The result is that the calculated efficiency of back shield- 
ing is underestimated, since the value of the experimental 
susceptibility is determined by applying finite increments 
AH. 



Figure 6: xfc calculated with Eq. (II. 6) from the values 
of Xg and Xj(0) (see text). 



II-B Sample B 

Sample B was machined from one of the original cylin- 
ders, in form of a parallelepiped of dimensions approx- 
imately 3x3x6 mm^ . Its calculated volume is « 
52.6 mm'^ and its demagnetizing field coefficient for a 
longitudinal field is N Ki 0.19. In a longitudinal field, 
its calculated moment for perfect flux expulsion is given 
by = 5.1 ± 0.2 • 10-3 X H cw?-G 

Measurements of the initial magnetization at 10 K are 
in fair agreement with this value. For Hdc above 3 G 
and up to 30 G AM / AH reaches a stable level about 
3.210 



-3 , 



which corresponds to the response of the 
grains alone. With the density ratio of 88% for this sam- 
ple and using Eq. ( |II.2| ) one finds / = 0.46 , yielding 
A = 0.19 r , i.e. the same value as derived for sample A. 

The shielding susceptibility was measured in this sam- 
ple by using a more sophisticated method, in order to 
reduce the effect of non linearity. After cooling the sam- 
ple at the working temperature in the d.c. field, the field 
was increased by 5 successive steps AH , and AA^ was 
measured. At the lowest fields, AH — 10 mG and (to 
keep a good signal/noise ratio) AH = Hdc/50 at the 
highest ones. Then, the value of AAl„/^,jAiJ was 
extrapolated to AH = by least square fit. 

Like in the case of sample A, all curves at Hdc < 10 G 
merge at high temperatures to a common curve which 
corresponds to the flux expulsion by the grains. The 
main difference with the sample A is that in the sam- 
ple B the onset of Josephson currents shielding occurs 
at higher temperatures. This is consistent with the fact 
that sample B is more dense, resulting in a better cou- 
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pling between grains; moreover, its size is larger, which 
also increases the total shielding magnetization. At low 
temperature, the magnetization curve at Hdc — 2QG 
reaches a level slightly above 60%, which corresponds 
to the low temperature level for the grains. 

The shielding response of the Josephson currents is 
obtained with the procedure already used for the sample 
A. Here the demagnetizing factor cannot be neglected 
{N K, 0.19). Two kind of quantities are to be consid- 
ered: (i) the responses Xg ^^id x of equivalent sam- 
ple without demagnetizing field (e.g. an infinitely long 
cylinder with the same cross-section for instance); here 
Xg is the response of the system of grains alone, without 
intergrain currents, and x is the total response of the sys- 
tem of intragrain plus intergrain currents, and (ii), the 
measured responses Xg and x ; they correspond to the 
measured moment for each case, normalized to the mo- 
ment for total flux expulsion in the volume of the sample. 
The relation between both kinds of quantities is given by 
A relation of the same kind holds 



M _ X _ X 
H ~ l-N ^ l+4TrNx' 

for Xg and Xg- With the use of Eq. (114), we finally 
obtain: 



(II.7) 



where Jt = l+Airx , Jig = l+'iirXg ■ Similar to the case of 
the sample A, an approximate curve has been determined 
for Xg by interpolating between the small Hdc curves at 
high temperatures, and the curve at Hdc = 20 G at low 
temperatures. Then the values of Xj have been derived 
from Eq. (II. 7) and plotted on Fig. 0. The set of curves 
is similar to the set for sample A, except for the higher 
onset temperature of the intergrains currents. 



Figure 7: Josephson currents shieldin g sus ceptibility as 
derived from the data and use of Eq. (|ll. 



Field Cooled magnetization data, normalized to the 
value of the moment for full flux expulsion, are reported 
in Fig. as a function of fields up to 30 G . 



Figure 8: F.C. (Meissner) susceptibility of sample B nor- 
malized to the moment for total flux expulsion, as a func- 
tion of field. 

At the lowest field and temperature, the F.C. mag- 
netization does not exceed 28% of its value for full flux 
expulsion. Furthermore, at low temperatures the curves 
versus field present a second maximum at about 5 G . We 
expect that this complicated behavior is due to the back 
shielding effect of the intergrain currents, as discussed for 
sample A. To take them into account, a relation similar 
to Eq. (II. 6) ( but with the demagnetizing effect taken 
into account) should be derived. The internal field is 
given as usual by Hi = H — AttNM , and the value of 
the local field seen by the currents is Hi = Hg Hi . Thus, 
under a temperature variation dT , 

^ ^ ^ (i^ _ 4.7VAf ) 



47riV fig 



dM 

It 



dT dT 

With dB/dT = jij dHi/dT , and using the relation 



M 



XfcH 
1 -N 



1 

An 



It, 



' d{B - Hi 
dT 



dT 



one obtains after integration: 



Xfc 



1 - N 
Attn 



(l-exp(-47r7V/)); 



1-N [l-iig 
(II.8) 

Here /ij = 1 + Anxj , with Xj reported on Fig. whereas 
the value of /ig was obtained using the relation jig = 
(1 — N)'jlg/{1 — Njlg) from the value of Xg as derived 
above. 
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The values of Xfc for _ff — > have been calculated 
using the values of Xg as determined above, and the val- 
ues of Xj &t Hdc = 0. The calculated value of Xfc 
found to be about —0.35 at T = lOK, whereas its mea- 
sured value was about —0.28. The discrepancy between 
measured and calculated values is larger here than in 
corresponding results for sample A. We believe that the 
origin of this discrepancy is the same as in the case of 
sample A, i.e. it stems from the nonlinear response ef- 
fect. This effect is numerically larger in sample B since 
here the onset of Josephson currents occurs in a range 
of temperature where Xg still varies strongly, contrary to 
the case of sample A. 

The above analysis shows (irrespectively to the above- 
mentioned discrepancy) that the back shielding effect 
leads to a strong reduction of the field cooled suscep- 
tibility as compared with the susceptibility of the grains 
alone. It is then easy to understand the complex behav- 
ior of Xfc ^s a funtion of field: at 10 K for instance, 
the onset of back shielding occurs at about 20 G , and 
its amplitude increases with decreasing field due to the 
increase of Xj ■ Starting from the two-step behavior of 
Xg expected from the theory of Wohllebeen et al. [291] 
(and seen in the data of sample A, where back shielding 
is less important), back shielding results on the double 
maximum shape of the measured curves. 

Ill Detailed Study of the Joseph- 
son Network Response. 

III-A Determination of the global criti- 
cal current. 

In this subsection we will present the procedure we used 
to extract the value of the macroscopic critical current in 
our sample B. This procedure is not quite trivial since we 
are interested in the dependence of the critical current 
on the background d.c. field in the sample, so we need 
to analyse the magnetization curves which depend both 
on the cooling field Hd.c. and on the field variation SH. 

Magnetization has been recorded at 10 and 20 K as 
a function of increasing AH with the smallest possible 
field steps {SH = lOmG), and starting from several F.C. 
states. From the AA4 data, it is possible to derive the 
value of the current response AAij as a function of AH. 
For this, we use Eq. (|ll.7l ) which can be written as : 



AM, 



AM - AMg 
{l-N-p)ji„ 



(III.9) 



grains system response is approximately derived in the 
same section: AMg w 3.2 10~^ x H cm^ — G at 10 X and 
AMg « 2.9 10^3 xiJcTO^-G at 20 if. Calculated values 
of AMj at 10 K are plotted in Fig. ^. Note the analogy 
of these results to the magnetization curves of classical 
type II superconductors with strong pinning (the differ- 
ence is that here AH plays the role of H). 



Figure 9: Shielding moment of the Josephson currents 
after cooling the sample at 10 X in a d.c. field in the 
range 0- 10 G. 

After cooling the sample at zero d.c. field, its response 
is obviously symmetric with respect to AH. When it 
is cooled in a finite d.c. field, it is not the case any- 
more, as was explained in the previous section. The 
magnetic moment of the sample just after cooling is 
Mfc ~ Mg + Mf-^ where Mf-^ is the positive mo- 
ment due to the back shielding by the Josephson currents 
which ha ve b een developed during the cooling process 
(see Eq. (IL4 and ILV)). So, the total moment produced 
by the intergrain currents after increasing the field by 
AH is Mj = M^^ + AM J. It is this moment which 
vanishes when Jc ^ (at large enough AH), and thus 
AMj approaches —Mf^. In Fig. |l^ we show the data 
recorded a.t T = 10 K and Hdc = 2 G. Curves recorded 
at positive and negative AH both converge to the value 



corresponding to --M^ : at 10 K 



where AMg is the magnetization of the grains alone; 
JI and Jig are defined in Section II-B. The value of the 



-Mf-^ is about 
1.1 10~^eTOu. 

When AH > Hdc , it is natural to expect that the 
response of the Josephson network does not depend on 
the initial state. A simple illustration can be given by 
analogy with Bean-like pinning in type II superconduc- 
tors |l^ . At large AH , when the induction profile has 
penetrated up to the center of the sample, the magneti- 
zation no longer depends on AH but only on Jc. If, as it 
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approximate evaluation is no longer relevant. 



Figure 10: Shielding moment of the Josephson currents 
after cooling the sample at 10 if in a d.c. field Hdc = 
2 G . Data are for positive and negative field steps. 



is the case in real materials, Jc varies with the induction 
in the sample, the magnetization depends on the total 
H , whatever the value of Hj^^ in which the sample was 
cooled. Actually, when plotted as a function of the total 
field Hdc + , the curves giving the total moment of 
network currents Mj'-^ + AA^^ merge in their "large" 
field part (i.e. above their maximum). The values 
have been calculated, with —M.J^ = 1.1 10~'^emM and 
0.710~'^eTOu for Hdc = 2G and 4G respectively. In or- 
der to obtain an optimal overlap between the curves, the 



following values have been used for ^M.g '■ 3.25 10 x 
Hcm^ - G at Hdc ^ OG, 3.22 lO'^ x H cm^ - G at 
Hdc = 2G and 4G. Indeed, the calculated values for 
AA4j at large AH are extremely sensitive to those for 



AM, 



This allows us to refine the determination of 



AMg . Note that the values quoted above do not differ 
by more than 1%, which is compatible with experimental 
accuracy and the possible variations of grains response 
with Hdc ■ 

Finally, from the knowledge of the true Josephson 
shielding response, in "large" fields we can now derive a 
rough evaluation of the critical current. Namely, above 
the maximum of AA4j , we calculate the value Jc of the 
average critical current which would give the value of 
the measured moment by use of the Bean formula |ig|] 
in a cylindrical geometry. For strong penetration, the 
magnetization is given in e.m.u. by M = JcR/3. With 
R = 0.15 cm and the values of the moment measured 
at 10 K and 20 K with Hdc — G , we obtain the data 
Note that the data are limited to 
H* — AnJrR below which the above 



displayed in Fig. |ll 
fields such that H « 



Figure 11: Calculated values of the averaged critical cur- 
rent Jc as a function of total field for strong field pene- 
tration. The big square corresponds to the initial Jc as 
determined in III-B 



III-B Low Field d.c. Response. 

We can now concentrate on the behavior of the Joseph- 
son currents moment at small AH . For this discussion, 
the currents susceptibility AMj/AH is plotted versus 
AH at 10 K and 20 K in Fig. |l| and |l|, respectively. 

At 10 i^T , after zero field cooling or cooling in a small 
field Hdc — 0.1 G , the response varies linearly with AH 
for small values of AH up to about 0.5 G . This lin- 
ear slope of AMj/AH is considered as the result of 
classical Bean-like pinning with critical current density 
Jc = H* /AttR , where 1 /AnH* is the initial slope of the 
curve |19|| . This initial slope is reported on the figure as 
the short dashed line which corresponds to H* = 2 G, 
leading to Jc « S.TA/cm^. 

At larger AH , the behaviour of currents suscepti- 
bility AMj/AH deviates from linear, which is the re- 
sult of both the magnetic-field dependence of the crit- 
ical current Jc , an intrinsic effect, and of the increas- 
ing degree of flux penetration into the sample , a purely 
size-dependent effect. Usually one uses the Bean model 
(generally with some i?-dependent critical current) in 
an appropriate geometry in order to deconvolute these 
two effects. However, one should keep in mind that 
the Bean model is a severe simplification of the prob- 
lem of constant pinning force, corresponding to the limit 
A ^ (i.e. the London penetration depth is supposed 
to be negligible with respect to the Bean penetration 
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Figure 12: Josephson currents susceptibility at 10 K vs. 
applied variation AH of field, after cooling in d.c. field 
Hdc ■ The meaning of dashed and dot dashed lines is 
explained in the text. 



Figure 13: Josephson currents susceptibility at 20 K vs. 
applied variation AH of field, after cooling in d.c. field 
Hdc ■ 



length). For the simplest sample shapes (thin slab or 
cylinder) it means that the condition A <C i? should be 
fulfilled, which is usually the case. However the situ- 
ation is more complicated for samples of square cross- 
section (like our one), where the effect of corners may 
become important even at A ^ i?. For such a geom- 
etry, the use of Bean model leads to exactly the same 
relation between critical current, external field and mea- 
sured magnetization as for the cylindrical ones, whereas 
one expects some difference if finite-A corrections are 
taken into account. At the present stage, we are not able 
to evaluate these corrections and therefore the values of 
the magnetization corresponding to our experimental ge- 
ometry with non negligible A . Nevertheless, we expect 
that it lies between the curves for two extreme limits. 
The upper one corresponds to the A ^ limit, where 
the magnetization is given simply by the Bean's formula 
for the cylinder: AttM/H = -1 + H/H* - H^/3H*^ 
for H < H* and AttM/H = -H*/3H for H > H* . 
A lower limit (thought rather artificial) consists of the 
"double slab" case in which the variation of magneti- 
zation is counted twice (once for each pair of edges): 
AttM/H = -1 + H/H* for H < H* /2 and AnM/H = 
—H*JAH for H > H* /2. Both curves are plotted in the 
Fig. |2| (dot dashed and long dashed curves respectively) 
for Jc = 3.7A/cm'^ and AMj/AH = -5.05 lO^^^m^ at 
AH 0. 

Let us now discuss the data starting from those ob- 
tained for low d.c. fields, Hd.c. = 0, O.IG. One can see 
that, after the initial linear part, the absolute value of 
the measured susceptibility is always smaller than the 
calculated ones. This corresponds to the decrease of Jc 
with increasing induction, as it is classically expected 
in granular materials, due to the suppression of inter- 
grain critical currents by magnetic field penetration into 
the Josephson junctions|^. This "classical" behavior for 
granular superconductors is usually analyzed by consid- 
ering the volume-averaged Josephson medium as a kind 
of type II superconductor in the dirty limit, provided its 
macroscopic penetration depth Aj is large as compared 
with the grains size 0, 

At Hdc > 2G the behavior of AMj/AH is quite dif- 
ferent: there is no initial linear slope, but a monotonic 
curvature is present down to the smallest AH . It is 
no longer possible to adjust a Bean like curve to the 
data. For instance, the Bean curve plotted on the lowest 
AH data for Hdc = 2 G is reported on the Fig. [l^ as a 
dashed line. It corresponds to a very small critical cur- 
rent of order 0.2 A/ cm? , and it is evident that the effec- 
tive screening current becomes much larger with increas- 
ing AH . Here, contrary to the case of Hdc = OG , the 
absolute value of the measured susceptibility is always 
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larger than the calculated one for a constant shielding 
current corresponding to the limit /S.H — > . This means 
that, whereas at Hdc = the effective screening current 
density stays constant and then slowly decreases with in- 
creasing AT? (which corresponds to classical Josephson 
pinning), at Hdc >2G it increases with A_ff suhlinearly 
(since a linear increase would correspond to a susceptibil- 
ity independent of AiJ ) . Such behaviour is quite unusual 
within the commonly accepted picture of screening in su- 
perconductors; indeed, we know that, for vanishing field 
excitations, the screening current may be either i) linear 
in Aff and reversible, as in the London (or Campbell 
pO|) shielding regime, or ii) constant (equal to the initial 
critical current Jc) and irreversible as in the case of the 
Bean-type critical state (or of any other known critical 
model, e.g Kim model exponential model |Q, etc). 

The above anomalous screening behaviour is even 
more pronounced at 2QK where, even after zero field 
cooling, no initial linear slope of AAij/AH can be seen 
in the data. All curves show the same anomalous be- 
havior as the data at 10 K in fields from 2 G . This spe- 
cific behavior is emphasized by plotting the difference 
between the measured susceptibility AAij/AH and its 
value for total flux expulsion AMj{0)/AH, versus AH 
on a Log-Log scale. In such a plot, at least in the regime 
of weak penetration, i.e. where AA4j/AH is larger than 
0.8 AAij{0)/ AH , sublinear variation of the shielding 
current density results in a logarithmic slope smaller 
than 1 for the curves of AM-j/AH (for AMj/AH 
smaller than 0.8 AA^j(0)/A_ff , we are in a regime of 
strong penetration where it is no longer possible to relate 
simply the variations of the moment response to those of 
the shielding current). In Fig. we have reported the 
three curves for which data are found in the range above 
0.8AMj{0)/AH, i.e. at T = 10K,Hdc = OG and 2G, 
and T = 20KHdc = OG. 

At 10 K and Hdc — OG , the logarithmic slope is about 
1 as expected, although at the smallest fields the curve 
crosses over to a smaller logarithmic slope closer to 0.5. 
At 20 and Hdc — OG the logarithmic slope is about 

0. 4 at the lowest AH . Approximately the same value 
of the slope characterizes the data obtained at 10 K and 
Hdc = 2 G , although the dispersion of data points at 
lowest AH makes its accurate determination difficult. 

The above anomalous behaviour makes it tempting to 
try a simple Ansatz for the behavior of the response cur- 
rent density of the system versus induction variations. 
Let us suppose that J esc AB" with a between and 

1. The case with a — 1 corresponds simply to classi- 
cal screening with penetration length A (since J c>c AB); 
the case with a = corresponds to constant J, i.e. the 
classical Bean case. Anomalous response arises for non 



Figure 14: Difference between the measured susceptibil- 
ity and its value for perfect shielding for selected data 
at 10 K and 20 K . The short dashed line represents a 
logarithmic slope 1 expected for a Bean critical state. 



integer a . For very small excitation AH , the length 
of induction penetration is small as compared with the 
size of the sample and we need to consider the effect of 
the excitation in the lowest order in AB only. For the 
purpose of illustration wc consider the simplest slab ge- 
ometry. Then the induction profile is determined by the 
Maxwell equation 



dB 
dx 



( AB 



\ABi 



(III.IO) 



where x is the coordinate perpendicular to the edge of 
the sample. For an external field AH , the induction in 
the sample is given by 



AB{x) 



(1 - Q!)47rJi 



ABi'- 



1 / 

- (XH 



1/(1-") 



(III.ll) 



where xh is the coordinate of penetration and Ji and 
ABi are normalizing factors; for a; = , AB = AH, i.e. 
XH = ABi°/(47r Ji(l - a)) ■ AH^'"' . Then, integrating 
the field profile (III.ll) over x, we get 



47rAM -I- AH 
A^^ 



AH^- 



(III. 12) 



where AM = AAi/V is the mean magnetization varia- 
tion due to the field variation AH. 

If we now compare the result ( III.12| ) with the data 
shown in Fig. 14, we find values of a in the range 0.4-;-0.5 
at both 10 K and 20 K. 
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Thus a simple choice for the relation between the 
screening current Jc and the induction variation Ai? al- 
lows us to imitate the experimental results for the sim- 
plest protocol of a weak monotonic AJ7 variation on top 
of a homogeneous state of the network. Nevertheless, it 
is evident that A_B has no clear meaning if the variation 
of H is non-monotonic or if the initial state is obtained 
by non-zero field cooling. Indeed, in the later case, in- 
duction in the sample varies during cooling due to the 
variation of /i^ with T , giving the response AAl j^"-^ as 
seen before. Furthermore, we will see below that the 
response is irreversible even for extremely low exitation 
fields. 



plays a log-log plot of the out-of-phase susceptibility x" 
versus the amplitude of the a.c. field, and for several val- 
ues of the F.C. static field. The response is irreversible 
down to the lowest a.c amplitudes, and the irreversibility 
increases with the superimposed d.c. field. All curves 



III-C Irreversibility: Very Low Field, 
low Frequency a.c. Response. 

Problems of sensitivity limit the range of small excita- 
tions which can be used in d.c. experiments. The pre- 
ceding results clearly show the sublinear nature of the 
low field response, but they do not allow its precise de- 
termination. In order to extend by several orders of mag- 
nitude the range of our lower excitations investigation, 
we have been led to perform a.c. susceptibility measure- 
ments. The use of a.c. response measurements is always 
questionnable when equilibrium (or quasi-equilibrium) 
properties are under investigation, since the results can 
be affected by the time-dependent part of the response 
function. It has been shown that the latter is the re- 
sponse of a very good conductor with complex conduc- 
tivity |36, Hence, it is necessary to work at low 
frequency, in a range where the susceptibility is roughly 
frequency independent. 

We present here preliminary results obtained on a long 
cylinder obtained by stacking several of the original sam- 
ple B cylinders. Measurements were done at 20 at a 
working frequency of 1.7 Hz in the equipment used for 
noise experiments The sample was simply shifted 

into the upper half of the third order gradiometer. At 
this temperature and frequency, we have verified that the 
in-phase susceptibility is almost frequency independent, 
which ensures that the results are mainly dependent on 
the (quasi) static part of the response. The suscepti- 
bility was recorded using classical method of SQUID 
magnetometry. We used a.c. excitation fields in the 
range 3 10~^-30 mG and the sample was cooled in d.c. 
fields from to 8.8 G. From the data, the values of the 
Josephson network susceptibility was extacted using the 
method developed in Section II, with the susceptibilities 
in Eq. ([L7) being complex quantities. The susceptibil- 
ity measured at A.2 K at the lowest a.c. amplitude was 
taken as the level for perfect diamagnetism. Fig. 15 dis- 



Figure 15: Out-of-phase susceptibility at 1.7 Hz as a 
function of a.c. field amplitude. 



Figure 16: In-phase susceptibility at 1.7 Hz as a function 
of the power of a.c. field amplitude H^-^ . 

follow a power law , with the same exponent close to 
0.5 . Going towards the smallest a.c. excitations, they 
show some downward bend which could be related with 
the approach to a linear regime (with Xj — 0) below 
0.1 toG, although the dispersion of the data is too high 
to conclude. The in-phase susceptibility x'j is plotted as 
a function of H^^,? in Fig. Here as well, the anoma- 
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Figure 17: Plot of 47rx" as a function of 1 + 47rXj- for the 
different values of Hdc ■ 



lous nature of the response is clearly seen. ^T^x'j behaves 
like (— 1 + (5 + '^H^'c) where both the constant 6 and the 
slope 7 increase with increasing superimposed static field 
Hdc ■ The dependence of the harmonic susceptibility on 
the a.c. field amplitude is a genuine proof of the existence 
of static irreversibility in the response. This is not as- 
tonishing by itself, but these results stress the anomalous 
aspect of this irreversibility. For instance, in the classi- 
cal Bean case with a weak penetration, it is known that 
1 -I- 47rX j and x'j are proportional to Hac whereas Fig. |l^ 
and |l^ clearly show the proportionality to Hl^^. A fur- 
ther evidence is provided by plotting Xj versus 1 -I- 47rXj 
as displayed in figure |l^. It can be shown that if the a.c. 
response is driven only by static irreversibility, both are 
proportional. In the Bean case, the coefficient of propor- 
tionality is 4/37r . In the figure ^ the part of data which 
lie in the range of 20% variation of x'j (where the rela- 
tions for slab geometry are approximately valid) show 
that x'j is indeed proportional to 1 -I- 47rXj, but with a 
bit smaller coefficient « 0.28 ± 0.03 . 

In order to understand the meaning of the above re- 
sults, we generalize the crude ad hoc model of Section 
III-B to the irreversible case. In order to do it, we gen- 
eralize the protocol of the Bean model. Namely, in the 
Bean model, the current is given by a step function of 
the variation of induction, J — JcSgn{AB) according 
to the sign of AB, as long as the induction variation is 
monotonous. If the sign of variation of B is reversed, 
J also changes sign, which can be written in terms of 
the variation of the current density (with respect to the 
initial current distribution obtained after monotonous 
variation of the field, Jinit , AJ = -2JinitQ{- ABnew) 



where Q{x) = {l/2){l + sgn{x)) and ABnew = B-Binit ■ 
Such a representation (which is not needed in the analy- 
sis of the Bean model itself) will allow us to construct the 
necessary generalization of the relation between current 
and variation of the field used in Eq.( III.IC). Actu- 
ally our goal here is rather limited: we are going to find 
a consistent description of the simplest hysteresis cycle 
which consists of the initial increase of AB from zero to 
ABinit , then reversing the sign of the field variation until 
the value of AB = —ABinit is reached, and then revers- 
ing dB / dt once more and finishing at AB final = ABinit- 
The description of this cycle will be consistent if we find 
that the value of the current density at the end-point, 
J final, coincides with the one after the original increase 
of the field ABinit, Jinit ■ This simply means that the 
hysteresis loop is closed. It is easy to check that the 
above condition will be fulfilled by the following choice 
of the AJ [ABnew) dependence: 



AJ = 



I J \ril — a T / A -Brie 
-Sgn{Jinit)2 J I 



e(-ABn 



(III.13) 

where Ji and ABi have the same meaning as in 
Eq.( |[II.l"0 ). Actually the only difference between the 
Eq. ( 111.13 ) and the original used in the Eq.( |lII.ld ) is 
the coefficient 2^~". The Bean model limit then corre- 
sponds to a ^ 0, so the above coefficient approaches 2 
as it should. Then instead of Eq. ( III.1C| ) we obtain: 



dAB/dx = ±2^-°'AAB°' 



(III.14) 



where A = AnJi/ABf . The induction profile, induced 
magnetization and harmonic response are calculated in 
the Appendix. The main conclusions are that the fun- 
damental components 1 -I- 47rXj- and x'j are both propor- 
tional to ho^~°' , and that their ratio R = x'j /{'^ + 47rXj) 
decreases from 4/37r to when a goes from to 1. For 
a = 0.5 , we get (cf. Fig. ^ R « 0.25 , a value which is 
in good agreement with the data presented on Fig. |l^. 
Note that the degree of irreversibility (measured by this 
ratio) is similar (although a bit lower) to the one of the 
Bean model. It should be emphasized that the numerical 
coefficient in Eq. ( 111.13 ) was "fitted" in order to obtain 
consistent (i.e. closed) hysteresis loop; one can expect 
that an analogous equation describing current variation 
after some more complicated history of the field varia- 
tions will contain another (history-dependent) numerical 
coefficient instead of 2^~". 

It can be seen from Figs. |l^ and |l^ that A — 
AttJi/AB" increases with increasing ambient d.c. field. 
It is natural to expect a decrease of Ji with increasing 
Hdc ■ The increase of A with Hdc means that ABf de- 
creases more quickly than Ji when Hdc increases. The 
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presence and behavior of the constant 6 cannot be pre- 
dicted on the basis of the above simple model. In fact, 
the latter neglects the possibility of elastic displacement 
of flux lines under the action of the external applied field. 
Such an effect would result in a response analogous to 
the Campbell response due to the elastic displacement 
of vortices in their pinning potential in type II supercon- 
ductors |2^. Campbell's response is linear, so it would 
result in an a.c. field independent positive contribution 
to x', whose amplitude would be inversely proportional 
to the strength of the restoring force of the pinning po- 
tential wells. It is natural to expect that the pinning 
force decreases with increasing ambient d.c. field in our 
granular system, due to the reduction of the junctions 
critical currents. Hence, such an effect would give a pos- 
itive contribution to which would increase with in- 
creasing d.c. field. This corresponds rather well to the 
behavior of the offset S seen in the data. 



IV Comparison with an existing 
theory of gauge glass: frustra- 
tion at H = 

In this Section we compare the experimental results de- 
scribed above with the theoretical results available for 
the randomly frustrated Josephson networks. We start 
from a simple estimate for the mean energy Ej ~ Tilc/'ie 
using the experimental value of the low-temperature, 
Z.F.C. (T = IQK.Hdc = 0) critical current density 
Jc ~ i.! Aj cm? . Using the estimate oq ~ for 
the mean size of the grains, one could naively obtain 
/j, « JcOg ~ IfiA and the corresponding low-temperature 



Josephson energy iJ^"™^ 20K (this value was derived 
from Jc measured at T = lOK, but we do not expect 
much difference in the intrinsic Josephson energies at 
T = lOK and at T ^ since the bulk transition tem- 
perature in Lai.8Sro.2Cu04 is Tc ~ 32K). However 
such an estimate is in contradiction with the measured 
value of the glass transition temperature Tg « 29K. In- 
deed, let us assume that the mean coordination number 
(number of "interacting neigbours" ) Z in the ceramics 
is around 6, as for a simple cubic lattice. Then for the 
estimate of the relation between Ej and Tg one can use 
the simulation data j|, Q which give Tg « 0.5Ej{Tg) = 
0.5Ej^-{l — Tg/Tc), where we took into account the linear 
dependence of Ej onTc — T close to the bulk transition 
temperature. As a result, one gets 



i.e. a factor 30 larger than the naive estimate above. 
However we will show now that this discrepancy may be 
resolved if we assume that the current network produc- 
ing the measured critical current density Jc was actually 
strongly frustrated in spite of the absence of background 
d.c. field in this measurement. 

The macroscopic critical current density Jc for a 
strongly frustrated Josephson network was calculated 
m Ref. g within the mean-field approach (we are not 
aware of any calculations of this kind beyond the scope 
of the mean-field theory) . It was shown that frustration 
strongly reduces Jc as compared to its value Jq for an 
unfrustrated system, Jc/Jq = ^^(1 - TjTgfl"^, where 
the factor 7 k, 0.065 was obtained by numerical solu- 
tion of the slow cooling equations |3^, describing 
the evolution of the glassy state under slow variations of 
temperature and magnetic field. In the low-temperature 
limit, this relation amounts to a factor 25 reduction of 
the Jc value with respect to Jq. Correspondingly, the 
characteristic value of the critical current for an individ- 
ual junction will be obtained as Ic ~ 25JcCLq « 25^A and 
results in a Josephson couphng energy Ej^ « 500^, in a 



2Ta 



1 - Tg/Tc 



600K 



(IV. 15) 



fairly good agreement with the above estimate (IV. 15). 

The above estimates show that the network of Joseph- 
son junctions in Lai.8Sro.2Cu04 is frustrated even in the 
absence of an external magnetic field. A careful reader 
could question this conclusion since we have used some 
results from the mean-field theory which may be a poor 
approximation for a 3D gauge-glass. We believe, how- 
ever that the qualitative result of the above estimates is 
sufficiently robust because a strong reduction of Jc with 
respect to Jq should be a general feature of a glassy net- 
work, so that unaccuracy due to mean-field approxima- 
tion cannot compensate for a huge discrepancy obtained 
between £:^'»™'= and the estimate ( |IV.15D . Additional 
evidence in favor of the glassy nature of our system is 
provided by the similarity of the \ow-AH diamagnetic 
response at T = 20K with zero as well as non-zero Hdc 
, as described in Section III above, as well as the low- 
frequency noise data obtained in Ref. |Q on the same 
type of ceramics. 

What could be the origin of that frustration? We be- 
lieve that most probably it is the result of the d-wave 
nature of superconductivity in cuprates 1 15 and random- 
ness of the crystalline orientations in ceramics jl^, 
It was shown there that the form of the effective phase- 
dependent Hamiltonian for such ceramics is of the same 
form as in (1.1) except for the fact that the random 
phases at -B = are just or tt depending on the 
mutual orientation of grains i and j. Therefore such a 
system at i? = is equivalent to the XY spin glass, with 
the low-temperature state characterized by a completely 
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random orientation of phases 4'i, as in the gauge-glass 
model with uniformly random distribution of a^j's. 

Therefore the low-temperature state is characterized 
by the presence of randomly distributed intergrain cur- 
rents and, therefore, of the magnetic field generated by 
these currents. It means that the actual phases will 
contain contributions due to the self-induced magnetic 
field. Its relative importance is characterized by the ratio 
of the corresponding magnetic flux penetrating elemen- 
tary loops of the ceramics ^sf to the flux quantum $0, 
i.e. just by the parameter /3i = 2ttCIc/c^o where C is 
the characteristic inductance of an elementary loop [^. 
Estimating the elementary inductance as £ ~ 27rao/ig 
and using Eg. ( IV. l'^ ) to estimate Ic, we obtain 



/3i 



$2 



0.1 



(IV. 16) 



so the self-field (screening) effects are relatively weak, 
though perhaps not always negligible. 

It is also of interest to estimate the effective penetra- 
tion depth Acer of a very weak magnetic field perturba- 
tion SH into the ceramics. Roughly, the value of Acer can 
be estimated as clo/VI^l ~ 15/im. Another (hopefully 
more accurate) estimate can be obtained using mean- 
field results which allow one to express Acer via the 
critical current density Jc'. 



7 



1/2 



25/i77l 



(IV. 17) 



where we inserted (as compared with Ref. [g) /ig ~ 0.35 
and approximated the random nearest-neigbour network 
by a cubic lattice with coordination number Z = 6, 
which amounts to the relation = Oq/Q between the 
effective interaction range ^0 and the intergrain distance 
ao- 

The characteristic magnetic field variation producing 
the critical current density Jc at the boundary can be 
estimated as AHc ~ AnXcerJcl c ~ 15mG, whereas the 
numerical solution [|j gives 



AHr 



2c7 



47rAcer>/c ~ SOmG 



(IV.18) 



Within the theoretical approach of Ref. Q , AHc marks a 
crossover between reversible (although still non-linear at 
AH < AHc) and irreversible penetration of the mag- 
netic field into the intergrain network. The value of 
AHc obtained in Eq.( IV.l"8| ) is on the lower border of 
the range of the field variations used to measure our d.c. 
magnetization curves, so we could just conclude that we 
always have SH ^ AHc and thus are producing the 
Bean-like critical state. Indeed, the data at Hdc = 0, 



T — lOK look compatible with such an interpretation 
(cf. Fig.([l^), where some deviations from the logarith- 
mic slope 1 (which is the characteristic of a Bean state) 
are seen at lowest AH < 50mG). 

However, as far as the data obtained at 10 K with 
d.c. fields Hdc — 2G , or all data at higher temperature 
(T = 20K), including d.c. and a.c. results at zeTO-Hdc, 
are concerned (cf. Figs.( p^| - |l^ )), the low-field magneti- 
zation response is drastically different from Bean-type 
predictions, as explained at the end of Section III. Qual- 
itatively, the most surprising feature of these data is the 
existence of a very broad range of AH within which the 
response is non-linear but still not like the critical-state 
one. We are not aware of any microscopic theory which 
predicts fractional-power behaviour of the shielding sus- 
ceptibility over such a broad range of AH variations. 
It cannot be excluded a priori that such a behaviour is 
related to a very wide range of intergrain critical cur- 
rents, which might exist in ceramics (till now we have 
neglected inhomogenity of intergrain coupling strengths 
in our theoretical discussion). Moreover, we may expect 
that the relative importance of such inhomogenities in- 
creases with the field and/or temperature (cf. Ref. |^ ) 

In Section V, we will try to formulate a new phe- 
nomenological model appropriate for the understanding 
of our data (leaving its theoretical justification for a fu- 
ture study); this model will be seen to be an interpolation 
between Campbell's and Bean's regimes of flux penetra- 
tion into hard superconductors. 

V Fractal model of diamagnetic 
response 

We showed at t he end of Section III-C that a simple gen- 
eralization, Eq.( III.14 ), of Bean's relation between vari- 
ation of the applied magnetic induction AB and current 
J results in reasonably good agreement with our data. 
However, contrary to the original Bean relation, the new 
one was not based on any physical picture; it was just 
a convenient description of the data. In this Section 
we propose a phenomenological model which provides a 
qualitative understanding of the irreversible diamagnetic 
behaviour mimicked by Eq.(III.14). 

We start from the picture of non-linear response of 
the current J to a variation of the vector potential 6A 
derived in Ref. |^ within the mean-field approximation 
and presented in Fig. 2 of that paper. Here the current 
induced by a variation of SA is linear at very small (5A, 
then grows sublinearly, and finally reaches its maximum 
value Jc at the critical 6 Ac such that the differential re- 
sponse {dJ/dA)sAa 0. At 6 A > SAc the numerical in- 
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stability of the slow cooling equations was detected and 
interpreted as an indication of the absence of any so- 
lution which would interpolate smoothly between zero 
and large (i.e. ^ SAc) values of 6A. In other terms, 
some kind of "phase slip" was expected to happen in 
the model ||], leading to a new metastable state, which 
would have lower (free) energy at the new value of the 
vector potential A' = A + 5 A (in other terms, a state 
similar to the one obtained by the F.C. procedure at 
constant A', which does not carry macroscopic current). 
Further increase of 6 A' — A~A' again induces a macro- 
scopic current until it reaches the maximum value Jc at 
SA' = SAc, and so on. Thus the whole J {6 A) dependence 
emerging from the mean-field solution Q is periodic; it 
leads immediately to the irreversibility of the response, 
since the inverse function dA{J) is multivalued: differ- 
ent vector potential values may correspond to the same 
value of current. Of course, such a periodic J{SA) de- 
pendence does not correspond to the usual Campbell - 
Bean picture, which would better be represented by 

Jc.-B.{SA) = J{6A)d{6Ac-6A)+Jc9{6A-SAc) (V.19) 

It is important to note that the J{6A) dependence 
Ref. Q was obtained from the space-independent so- 
lution for the glassy correlation function Qjj{t,t') —< 
cos{(j)j{t) — 0j(^')) >i such an approximation, being 
reasonable for the description of smooth "adiabatic" 
transformations in a system with long-range interac- 
tions, will probably break down when the jump from one 
metastable state to another happens. In other terms, the 
above-mentioned "phase slip" should have something to 
do with spatially inhomogeneous processes like vortex 
penetration in hard type-II superconductors. The prob- 
lem of the solution of the general history- and space- 
dependent system of integral equations (which may be 
derived following the method of Ref. Q) is formidable 
and the method to solve it is still unknown. Therefore we 
can only speculate on possible properties of its solution. 
The simplest idea would be that the macroscopic J{SA) 
response becomes (after averaging over inhomogenities of 
the space-dependent solutio n) sim ilar to the Campbell- 
Bean type of the response ( V.19| ). Indeed, our analysis 
of the low-field diamagnetic response at T = IQK and 
Hext = (Section III-B) developed in Section IV on the 
basis of such an assumption, is in reasonable agreement 
with the data. However other sets of data (for higher 
temperature and/or lower field) are described by com- 
pletely different Ansatz ( III.13| ). We will now propose a 
(phenomenological) generalization of the J{6A) relation 
compatible with Eq.([II.13). The relation we are looking 
for should be an intrinsic (i.e. independent on the sample 
geometry) and general (i.e. usable for an arbitrary mag- 



netic history of the sample) relation between the current 
and variation of the vector potential. Remember that 
Eq. (III. 13) was written for the simplest nonmonotonic 
variation of AB, and that it relates the true vector J 
and the pseudovector AB. So, in writing this equation, 
some additional information on the geometry of the sam- 
ple has been used (we use the simplest slab geometry). 
Thus a natural basic equation should relate the current 
density J and the variation of the vector potential SA. 

In a generalized model, the diamagnetic current re- 
sponse should possess two major properties: 

i) it must scale as some fractional power a w 0.5 with 
the amplitude of exitation field SB, and 

ii) it must be strongly irreversible (as it follows from the 
analysis of the ratio 47rx"/(l - 47rx') « 0.28 ). We con- 
sider these two conditions in sequence. 

The condition i) is rather easy to fulfill: it is enough 
to suppose that the differential response of the current 
to the variation of the vector potential SA is given by a 
non-linear generalization of the London relation 



dA^ An^-ff^-^' 



(V.20) 



where the current-dependent "effective penetration 
depth" is given by 



(V.21) 



In the case of a monotonic field variation applied 
to an initially uniform induction distribution, the 
Eqs.( |V.20| , [V.21D lead to the simple relation J cx AB" 
with a = (l-fK)~^. Indeed, with dA — AB dx and 
approximating dAB/dx by AB/Xejf , one obtains J cx 
A_BT+« . Thus we need to choose k w 1 in order to re- 
produce the observed scaling with a w 0.5. 



However, the set of equations (V.2C, V.21) does not 
fulfill the second condition ii) above: the corresponding 
solutions are reversible, as it follows from the existence 
of a single-val ued function SA{J) (x which fol- 

lows from Eqs.( V.2C , V.2l[). I n other words, the system 
described by Eqs.( V.20 , V.21 ) would exhibit nonlinear- 
ity and harmonics generation, but would not show finite 
x"(w) in the uj ^ limit. In order to avoid this inconsis- 
tency, we need to formulate a model with the same kind 
of scaling between SA and J as in Eqs.( V.20 , V.21 ), but 
with a nonmonotonic J{SA) dependence allowing for the 
irreversible behaviour. 

A model obeying very similar properties was formu- 
lated and studied in Ref. |^ in a different physical con- 
text (one-dimensional spin glass). The low-energy spin 
configurations in this model are described by a phase 
variable f{x) G (— tt, tt) such that two such configura- 
tions (local energy minima) which differ by a phase shift 
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Sip{xo) = <f> in a region around some point a;o, have 
a characteristic energy difference E{^) cx <i)^/^ and a 
characteristic spatial extent of the phase deformation 
X($) oc This scahng holds for the intermediate 

range of phase deformations ipo <^ ^ <^ n; at smaller 
^ < '/'o ^ 1 the energy cost of deformation is oc 
whereas at $ ~ tt the energy growth obviously satu- 
rates due to 27r periodicity. The above £'($) scaling 
leads to a sublinear growth of the characteristic "force" 
/($) — dE/d^ cx with $ in the same intermediate 
range. The main contribution to the second derivative 
d^E/d^^ (curvature of the energy valleys) comes from 
the smallest scale $ ~ ipo, i.e. from the curvature of 
individual local minima. It was explained in Ref. |4l| ] 
that such a scaling means a fractal organization of the 
energy minima as a function of ip with fractal dimen- 
sionality Df — 1/3. It means that the number of energy 
minima discernable on a scale ip grows as A/" oc (p~^^^ 
at finer scales; new minima appear primarily due to the 
splitting the older (broader) ones. This picture emerged 
in Ref. |4l| from the microscopic analysis of the original 
Hamiltonian for a one-dimensional spin-glass model for- 
mulated in Ref. j42j . We can borrow the qualitative fea- 
tures of this construction for our present purpose (leaving 
for future studies the problem of its microscopic justifi- 
cation for the case of superconductive glasses). 

Suppose that the free energy F{SA) of the Josephson 
network behaves (as a function of vector potential varia- 
tions with respect to a "virgin" state with a homogeneous 
induction) in a way similar to £'(1^9) at ip n. Namely, 
suppose that the free energy is parabohc, 6F oc {SA)'^, 
at very small variations of vector potential SA < SAd, 
but on a larger scale, SA 3> SAd, it contains many lo- 
cal minima whose characteristic free energies scale (with 
respect to the lowest state with SA = 0) as 



F{SA) oc {SAf+^ forSAd < M < SA^ 



(V.22) 



with the exponent 9 E (0,1) (see the definition of SAc 
below) . Then the characteristic value of the current J = 
-dF/dA scales as 



Jchar{SA) ~ Jc 



( SA 



\SAd 



(V.23) 



in the same interval of SA. At large SA > SAc variations, 
the growth of the induced current should saturate at the 
true critical current value Jc, so we can estimate 



SAc^SAdiJc/JciY^'- 



(V.24) 



On the other hand, weak SA <^ SAd leads to the usual 
linear London (or Campbell) response with an effective 



penetration depth Ai; matching at SA ^ SAd leads to 
the following estimate: 



SAd ^ — Jci^i 
c 



(V.25) 



The estimate (V.23) looks very much like the previous 
version defined by (V.2C, V.21), so one can find the re- 
lation between the exponents: 



= 1/(1 + 2k) = a/{2 - a) « 0.3 



(V.26) 



Figure 18: Picture of a fractal J{SA) landscape. An 
example of a hysteresis loop is shown. 

However the whole picture is substantially altered: the 
current is now supposed to be an (irregularly) oscillat- 
ing function of SA (see Fig. |lj), thus only its envelope 
Jchar{SA) defined on a scale SA follows the scaling re- 
lation ( V.23| ). As a result, the inverse function SA{J) 
is multivalued and the irreversibility of the response is 
ensured. Similar to the spin-glass model of Ref. pl|] , 
the fractal dimensionality Df oi the low-energy valleys 
can be defined; it is given now hyDf = l — 9^^ 0.7. 
The proposed picture is based on the existence of two 
substantially different scales of currents, Jd and Jc, and 
corresponding vector potential variations SAd and SAc] 
thus it can be compared with the usual Campbell-Bean 
picture of critical currents in the same way as the thermo- 
dynamics of type-II superconductors is compared with 
that of the type-I ones. 

In order to describe quantitatively the diamagnetic re- 
sponse in the "fractal" range ( |V.22 ) we need to determine 
the distribution function V[ J{SA) ] (which would lead, in 
particular, to the estimate (V.23) for Jchar{SA)). More- 
over, in general, a relation of the type of ( V.23| ) could be 
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nonlocal (i.e. the current depends on the SA{x) distri- 
bution in some region of space, whose size may depend 
on 6 A itself (see again Ref. We leave this com- 

plicated problem for future studies, and just note here 
that merely the existence of relation (V.23) is sufficient 
for the existence of some "natural" properties of the re- 
sponse (like the presence of a closed hysteresis loop, as 
it was assumed in Section III-C). 



VI Summary and conclusions 

In this paper, we have presented experimental results on 
the low temperature {10 K and 20 K) response of the 
granular HTc superconductor Lai.8Sro.2Cu04 to small 
field excitations. The general properties of the magnetic 
response were investigated in two samples (A and B) 
differing by the strength of the coupling between grains. 
By cooling the samples in various d.c. fields up to 20 G 
and applying small field increases, we were able to mea- 
sure the shielding response of the material and to derive 
a method, inspired by the work of Dersh and Blatter 
||2^ , to extract from the data the polarizability of the 
intergrain currents system. The field cooled (F.C.) mag- 
netization was measured in fields up to 20 G . Analysis 
of the results led to the conclusion that i) the structure 
of the grains is polycrystalline, resulting in a step de- 
crease of the F.C. magnetization with increasing field, 
which can be interpreted on the basis of the model by 
Wohllebeen et al. ii) self shielding (pinning) by 

the intergrain currents when lowering the temperature 
strongly reduces the value of the F.C. magnetization; iii) 
there is no macroscopic Meissner magnetization due to 
the system of intergrain currents. 

Further detailed study of the response of the Joseph- 
son network was performed in sample B. It was shown 
that the response is asymmetric with respect to the sign 
of variation of the applied field after field cooling; this 
is due to the shielding currents pinned during cooling. 
The macroscopic critical current is found to be strongly 
reduced by moderate values of the external d.c. field, 
about 2 G . 

Very low field magnetization measurements were per- 
formed by applying field steps of 10 mG or low frequency 
a.c. fields in the range 50 fiG to 30 mG, after cooling in 
d.c. fields up to 8.8 G. The results show that the re- 
sponse is strongly non linear, the shielding current grow- 
ing sublinearly with increasing applied field. Further- 
more, the a.c. results show that it is strongly irreversible 
down to the smallest excitations used. It is shown that 
a non-linear relation between the shielding current and 
the induction, J oc AB" with a ~ 0.5, together with a 



natural assumption about the existence of a closed hys- 
teresis loop, give predictions in a reasonable agreement 
with the data. 

Theoretical analysis of our experimental results was 
developed on the frame of the existing "gauge-glass" 
theories. It was show that the extremely low value of 
the low-temperature, zero-field critical current density 
(Jc « 3.7 A/ cm^ at 10 K) together with the rather high 
temperature of the transition to the low-temperature 
glassy state, can be coherently interpreted only under 
the assumption that the Josephson network is strongly 
frustrated even at zero applied field. This contradicts 
the usual assumption that frustration in the interactions 
arises only due to the local magnetic induction, but sup- 
ports the hypothesis of the existence of a large propor- 
tiont of TT-junctions in the granular system. These tt- 
junctions are possibly due to the d-wave nature of the 
pairing, combined with the randomness of grain orienta- 
tions in Lai.8Sro.2Cu04 ceramics. 

Finally, a new model of diamagnetic response in the 
glassy state of granular superconductors was developed 
in order to describe the anomalous (fractional-power) be- 
havior of the shielding current response. This model, 
based on the idea of a fractal organization of the free 
energy landscape in the granular network, can provide a 
qualitative account for the main features of the anoma- 
lous response. Its further development will be the subject 
of future studies. 
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A Appendix 

The hysteretic behavior of the current as a function of 
the induction variations is represented by the relation: 



A J = ±2 



(l-a)i/ 



Jr 



\AB\ 

Br 



(A.27) 



V — when starting from zero induction state, and 1 
otherwise. A J = J— Jo and AB = B — Bo where Jo and 
Bo are the (old) values just before the last reversal of the 
sign of variation of B . The Ansatz A.27 ensures that we 
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have a stable closed hysteresis loop, and that there is no 
hysteresis for a = 1 which describes the London case. 
Figure |l^ displays the shape of the hysteresis loops for 
three characteristic values of a . The induction profile is 



Figure 19: Hysteretic behavior of the current for a vari- 
ation of the induction between —Be and Be ■ 

determined by the Maxwell equation which leads, for the 
case of weak penetration, to 



dAB 

dx 



(A.28) 



where A = A.'kJcI B^ ; x is the distance from the edge of 
the sample. After increasing applied field from to /lo, 
starting from zero induction state, the induction profile 
is given by B~°'dB = —Adx , leading to: 

where 

B ^ {hi-" - {I - a)Ax)^ . (A.29) 

Field penetrates till x = Xh^ — — a)A. 

When h decreases from /lo, we get [Bo — B)-°'dB = 
— 2^~"dx . Hence: 



1 



2^-"A J. 



Ba-B 



r"d^ 



2i-"yl(l - a) 



((Bo - By-" - [K - hf 



Modification of induction relative to Bo extends up to 
Xh = {ho ~ hf-^ji'-^A . For < X < Xh , 



l-Q 



B = Bo-2 I ( ^^V-^ ) -{l-a)Ax\ (A.30) 



Figure 20: Induction profile as a function of applied field. 
Ro is the maximum abcissa of field penetration. 



where Bo is given by Eq. (A.29). When h 



reached, Eq. (A.3C) gives simply B — —Bo- After re- 



versing the sign of variation of h once more, the profiles 
are simply symmetrical of those given by Eq. ( |A.30 ). 
Examples of induction profiles are given on Fig. GOfor 
a = 0.5 . The average induction can be derived now. 
After some algebra, one obtains: 



<B>-- 



A 



2-a 
for ho > , and 

<B>= — 



X h 



2-a 



xhl 



a 



-1 



1 - h/ho 



1 + h/ho 



2-a 



(A.31) 



for ho <0 . 

For a sinusoidal excitation h = ho coswt, one gets 



<B> Ah\ 



ho 2 — a 

for 2n7r < ujt < {2n + l)n 
<B> Ahl-" \ ^ 

; = 2 

ho 2 — a 

for {2n — 1)tt < u)t < 2mT . 



1 - 2 



1 — cos ujt 



2-a 



1 + COS Ujt 



2-a 



(A.32) 



Since < B > /ho = 1+ < M > /ho, Fourier trans- 
formation gives the values of 1 + Atty' and 47rx" . This 
can be done numerically. Figure 21 displays the ratio 
47rx"/l + 47rx' as a function of a . 
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Figure 21: Values of 47rx"/l + ^ttx' as a function of a . 
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